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Abstract 
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'The  purpose  of  this  contract  is  the  development  of  a  Low  Noise,  Closed 
Cycle,  Joule  Thomson  refrigeration  system  for  10K  operation.  This 
report  summarizes  work  completed  to  date  on  the  detailed  thermodynamic 
analysis  and  physical  design  of  the  refrigerator.  A  complete  enthalpy 
balance  for  the  three  stage  device  is  presented.  Fabrication  of  the 
first  two  refrigeration  stages  (argon  and  hydrogen)  has  begun.  A 
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1.0  Introduction 


This  report  summarizes  the  results  of  an  engineering  design  program 
undertaken  by  MMR  Technologies,  Inc.,  for  the  Office  of  Naval 
Research  under  Contract  No.  N00014-85-C-0428.  The  purpose  of  this 
program  is  to  develop  a  low  noise,  closed  cycle,  Joul e-Thomson, 
microminiature  refrigeration  system  for  10K  operation. 

A  thorough  thermodynamic  analysis  has  been  completed  on  the  entire 
three  stage  refrigerator  and  is  summarized  in  Sections  2.1  and  2.2. 
Engineering  layout  of  the  gas  channels  has  been  completed  for  the 
first  five  layers,  the  argon  and  hydrogen  refrigeration  stages,  of 
the  proposed  eight  layer  refrigerator.  The  physical  design  and 
function  of  this  portion  of  the  cooler  is  described  in  Section  2.3. 

A  prototype,  single  gas,  two  stage  compressor  has  been  designed 
built  and  tested.  The  details  of  this  and  the  results  obtained  are 
described  in  Section  3. 

2.0  Refrigerator  Design 

During  the  present  reporting  period  (August  1  to  October  15,  1985) 
the  thermodynamic  design  of  the  ONR  helium  refrigerator  has  been 
completed.  The  results  of  the  detailed  enthalpy  balances  for  the 
argon,  hydrogen  and  helium  stages  are  described  in  Section  2.1.  A 
description  of  the  required  operating  conditions  of  the  heat 
exchangers  in  each  stage  is  included  in  Section  2.2.  Finally  the 
physical  description  of  the  argon  and  hydrogen  stages  is  given  in 
Section  2.3. 

2.1  Thermodynamic  Analysis 

The  thermodynamic  analysis  of  the  ONR  refrigerator  proceeds  from 
detailed  enthalpy  balances  of  each  stage  of  the  device.  These 
balances  are  solved,  one  after  the  other,  beginning  with  the  helium 
stage,  then  hydrogen  and  argon  taking  into  account  the  refrigeration 
capacity,  the  thermal  masses  cooled,  the  heat  leaks  present,  the 
enthalpy  differences  expected  in  the  gas  streams  and  the  theoretical 
heat  exchanger  efficiencies.  The  result  of  these  calculations  are 
the  gas  flowrates  required  in  each  stage,  and  are  summarized  in  the 
following  sections.  An  overall  enthalpy  balance  for  the  ONR  helium 
refrigerator  is  shown  in  Figure  2.1-1,  and  a  physical  schematic  of 
the  device  is  shown  in  Figure  2.1-2. 

2.1.1  Enthalpy  Balance  -  Helium  Stage 

The  enthalpy  balance  for  the  helium  stage  that  applies  during  the 
transient  cooldown  period  can  be  written; 

"He  “hg  *  I  Eqi  *  * 
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Figure  2.1-1 


where  the  denotes  the  thermal  masses  between  23K  and  10K,  qj 
denotes  the  heat  leaks,  and  t  is  the  length  of  the  transient  period. 
The  thermal  masses  and  power  dissipation  terms  are  tabulated  in 
Table  2.1-1.  The  thermal  mass  terms  for  sensor  and  refrigerator 
material  are  based  upon  the  physical  schematic  shown  in  Figure 
2.1-2.  The  heat  leak  is  that  caused  by  thermal  conduction  through 
the  refrigerator  material  (glass)  between  23K  and  10K.  The 
radiation  term  results  from  a  300K  ambient  exposed  to  the  gold 
flashed  surface  of  the  coldstage  (emissivity  assumed  to  be  0.02). 
When  the  enthalpy  balance  is  solved  the  dominant  term  for  the  helium 
stage  is  the  required  refrigeration  at  10K,  specified  to  be  10 
milliwatts.  The  volumetric  flowrate  of  helium  based  upon  at  rapid 
cooldown  of  three  minutes  is; 

VHe  =  1.3  liters/min.  at  STP. 


2.1.2  Enthalpy  Balances  -  Hydrogen  Stage 

stage  during  the 


(2) 


portions  of  the 
The  power 

dissipation  terms  (qj)  are  due  to  the  conduction  heat  leak  through 
the  refrigerator  material,  radiation  incident  on  the  hydrogen  stage 
and  the  precooling  required  by  the  helium  stream.  The  numerical 
values  of  these  terms  are  shown  in  Table  2.1-1.  For  this  stage 
there  is  no  single  dominant  term  that  determines  the  refrigerator 
performance.  The  refrigerator  thermal  mass  (qqiass)  and  the 
precooling  required  by  the  helium  stream  (qHe )  combine  to  account 
for  most  of  the  heat  load  in  this  portion  of  the  device.  The  result 
of  this  balance  is  a  volumetric  flowrate  of  hydrogen  of; 

V„  =  1.20  liters/min.  at  STP. 


The  enthalpy  balance  that  applies  to  the  hydrogen 
transient  cooldown  period  is  given  by; 

ViH1e  '  FEiq1  *  -qj 

£  l  J 

The  thermal  mass  terms  (q, )  are  taken  to  be  those 
refrigerator  cooled  by  hydrogen  from  9GK  to  23K. 


2.1.3.  Enthalpy  Balances  -  Argon  Stage 

The  enthalpy  balance  that  applies  to  the  argon  stage  during  the 
transient  cooldown  period  is  given  by; 

=  IEi  "f  +  Ej  ’j  <3> 

The  thermal  mass  terms  (qi )  represent  the  entire  material  of  the 
refrigerator  that  is  cooled  from  an  initial  ambient  temperature  of 
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300K  to  the  steady-state  condition  at  which  the  boilers  are  at  9QK 
and  a  gradient  exists  across  the  argon  precooler  from  90K  to  the 
ambient  of  300K.  The  power  dissipation  terms  ( q j )  are  the  result  of 
the  conduction  leak  through  the  refrigerator  material,  the  radiation 
incident  on  the  refrigerator  from  the  300K  ambient,  and  the 
precooling  required  by  the  hydrogen  and  helium  streams.  The  values 
of  these  terms  are  shown  in  Table  2.1-1.  As  would  be  expected  for 
this  high  temperature  stage,  the  thermal  mass  term  completely 
dominates  this  calculation.  This  is  due  to  the  relatively  large 
specific  heat  of  the  refrigerator  glass  at  high  ambient  temperature. 
The  result  of  this  balance  is  a  volumetric  flowrate  of  argon  of; 

V.  =  3.0  liters/min.  at  STP. 

Ar 


2.2  Heat  Exchanger  Operation/Estimates  of  Required  Approach 
Temperatures 

A  method  of  determining  the  feasibility  of  the  proposed  refrigerator 
design  is  to  calculate  the  average  approach  temperature  in  each  heat 
exchanger  section.  These  calculations  use  the  schematic  physical 
layout  shown  in  Figure  2.1-2  for  the  allowed  heat  exchanger  area, 
and  the  mass  flowrate  resulting  from  the  enthalpy  balances  of 
Section  2.1.  The  average  approach  temperatures  required  for 
steady-state  operation  (described  by  the  overall  enthalpy  balance) 
is  calculated  from; 


q  =  K  A  AT  =  mAH 
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Here  we  simply  say  that  the  precooling  achieved  in  a  given  section 
of  precooler  is  equal  to  the  product  of  the  thermal  conductivity  of 
the  exchanger  material  (K)  divided  by  the  exchanger  wall  thickness 
(d)  times  the  area  (A)  and  the  temperature  difference  (AT)  between 
the  two  gas  streams. 


So,  ATaverage 


mAHd 


(4) 


where ; 

m  =  mass  flowrate  through  the  exchanger 

AH  =  change  in  enthalpy  required  by  the  enthalpy  balance. 

A  =  area  for  heat  exchange 

d  =  thickness  of  exchanger  wall 

=  thermal  conductivity  of  exchanger  material  of 
construction 


K 


Precooling  Stage 


Helium 


Calculated  Maximum  Approach  Temperature 
Required  in  Precooler  Stages 
of  the  ONR  Refrigerator 


*  high  flow  calculation  based  upon  4000  psig  operation  for 
fast  cooldown. 

*  assuming  .003"  interlayer  design. 

*  approach  temperature  =  temperature  difference  between  two 
gas  streams  shown  in  narrative. 


Table  2.1-2 
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Table  2.2-1  shows  the  temperature  change  that  is  required  in  each 
heat  exchanger  section  of  the  ONR  refrigerator. 

For  example,  in  Table  2.2-1  we  see  that  the  argon  stage  must  precool 
the  hydrogen  stream  from  the  ambient  temperature  of  300K  to 
approximately  90K.  The  quantities  used  in  equation  (4)  above  are 
determined  as  follows; 

mH  ,  determined  from  the  enthalpy  balance  of  Section  2.1.2, 
n2  equals  2.0  x  10_,jgm/sec. 

AH  ,  determined  from  the  overall  enthalpy  balance  of  Figure 
2.1-1,  equals  2.88  x  10 3 joules/gm. 

d  ,  heat  exchanger  wall  thickness,  approximately  2.5  x  10~3cm 
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A  ,  heat  exchange  area  available,  approximately  4.8  cm 

K  ,  thermal  conductivity  of  the  heat  exchanger  material  (at 
approximately  15QK).  8  x  10"^watts/cm  K 

With  these  values,  equation  (4)  yields,  ATavera  e  =  3.1  K.  In  a 

similar  manner  the  required  approach  temperature  can  be  calculated 
for  each  section  of  the  cooler.  The  results  of  these  calculations 
are  shown  in  Table  2.2-2. 

2.3  Physical  Design/Function 

During  this  reporting  period,  the  design  of  the  first  five  layers  of 
the  ONR  helium  refrigerator  have  been  completed.  The  gas  channels 
have  been  layed-out  at  25x  true  scale,  rubilith  has  been  made  of 
these  patterns,  and  then  photo-reduced.  The  resulting  "true  size" 
layer  patterns  are  shown  in  Figure  2.3-1.  In  each  of  these  layers 
the  gas  channels  have  been  sized  to  accommodate  the  gas  flowrate 
calculated  in  Section  2.1,  and  to  produce  the  pressure  drop  required 
in  each  section  of  the  device.  The  theoretical  design  of  the  helium 
stage  has  been  completed  and  the  detailed  layout  of  these  last  three 
layers  will  begin  soon. 

The  operation  of  each  of  the  first  five  layers  can  best  be 
understood  with  reference  to  Figure  2.3-1.  Layers  1,  2,  and  3 
comprise  the  argon  stage.  As  shown  schematically  in  Figure  2.1-2 
the  gas  channels  of  this  stage  wrap  approximately  two-thirds  of  the 
way  around  the  outer  portion  of  the  1.5"  diameter  disk  of  the 
cooler.  High  pressure  argon  gas  enters  layer  2  through  a  small 
inlet  hole  located  at  approximately  8  o’clock.  The  warm  gas 
circulates  clockwise  in  this  layer  and  is  precooled  with  only  a 
slight  pressure  drop.  Upon  reaching  approximately  2  o'clock  in  this 
layer,  the  argon  undergoes  its  first  expansion  reaching  a  pressure 
of  approximately  20  atmospheres  at  point  "A".  At  this  point  the  gas 
stream  splits,  the  major  portion  passes  downwards  into  layer  1. 

Once  in  layer  1  this  stream  circulates  counter-clockwise  and  exits 
the  device  through  the  exhaust  hole  in  layer  1.  The  smaller 
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Gas  Stream 


Argon 


Hydrogen 


fraction  of  the  argon  stream  remains  in  layer  2,  undergoes  a  second 
expansion  and  liquifies  in  the  argon  boiler.  The  liquid/gas  mixture 
produced  in  the  boiler  of  layer  2  passes  upward  into  layer  3,  boils 
and  circulates  counter-clockwise,  finally  exiting  the  cooler  through 
the  exhaust  hole  in  layer  3.  By  the  arrangement  of  gas  channels  in 
adjacent  layers,  a  counter-current  heat  exchanger  is  produced.  The 
warm  stream  in  layer  2  is  "sandwiched"  between  two  cold  streams  in 
layers  1  and  3. 

Layers  4  and  5  comprise  the  hydrogen  stage  of  the  device.  Its 
operation  is  similar  to  that  of  the  argon  stage  just  described; 
however  the  gas  channels  are  sized  for  the  required  hydrogen 
flowrate  and  a  Hampson  cycle  is  used.  The  gas  channels  of  this 
stage  wrap  "once  around"  the  outside  of  the  refrigerator  disk.  As 
can  be  seen  clearly  from  Figure  2.3-1,  the  hydrogen  capillary  and 
boiler  extend  beyond  the  end  of  the  argon  stage  below;  thus 
providing  thermal  isolation  for  the  hydrogen  boiler  (at  23K)  from 
the  argon  boiler  (at  90K).  High  pressure  (60  atm)  hydrogen  gas 
enters  layer  4  through  a  small  hole  located  at  7  o'clock  in  that 
layer.  It  circulates  clockwise  in  the  precooler  channels  of  that 
layer  and  is  cooled  while  undergoing  little  pressure  drop.  After 
passing  over  the  argon  boiler  and  reaching  90K,  the  hydrogen  stream 
expands  and  liquifies  in  the  boiler  located  at  6  o'clock  in  layer  4. 
The  entire  gas/liquid  stream  passes  upward  into  layer  5.  Here  the 
cold  stream  passes  counter-clockwise  through  the  low  pressure  heat 
exchanger  and  finally  exits  the  device  through  the  large  hole  in 
layer  5. 

At  present  the  center  of  the  refrigerator  disk  is  not  used.  The 
helium  stage  will  occupy  all  the  area  of  the  disk.  The  final  heat 
exchanger,  expansion  capillary  and  boiler  ( 10K )  of  this  stage  will 
extend  beyond  the  hydrogen  boiler  (23K)  into  the  presently  unused 
center  portion  of  the  refrigerator  disk. 

3.0  Compressor 

A  prototype,  single  gas  (Air),  two  stage  compressor  has  been 
designed  ,  built  and  tested.  Its  design  was  based  on  the  design 
principles  and  mode  of  operation  outlined  in  the  first  progress 
report. 

An  electrically  driven  hydraulic  pump  delivers  hydraulic  fluid 
through  an  electronically  timed,  four-way  valve  to  a  double  acting, 
double  rod  hydraulic  cylinder  which  drives  the  first  and  second 
stage  pistons  of  the  pneumatic  circuit.  The  return  fluid,  after 
passing  back  through  the  four-way  valve  is  cooled  through  an 
air-cooled  radiator,  before  returning  to  the  reservoir.  The 
cylinder  bore  and  stroke  were  sized  to  give  a  cycle  speed  of 
fourteen  cycles  per  minute.  The  layout  of  the  system  is  shown 
schematically  in  Figure  3.0. 


3.1  Component  Description 

The  key  components  of  the  compressor  are  summarized  below: 

3.1.1  Hydraulic  Cylinder 

Double  acting,  double  rod  0-ring  sealed  piston  in  aluminum  housing. 

3.1.2  First  Stage  Cylinder 

Aluminum  cylinder  2.50"  diameter,  2.00"  stroke,  bearing  surface  hard 
anodized  and  honed  to  a  2-4  micro-inch  finish. 

3.1.3  Second  Stage  Cylinder 

Aluminum  cylinder  0.75"  diameter,  2.00"  stroke,  bearing  surface  hard 
anodized  and  honed  to  a  2-4  micro-inch  finish. 

3.1.4  Seals 

Single,  spring  loaded  lip  seal  of  filled  teflon  for  each  piston. 

3.1.5  Pistons 

Aluminum  pistons  anodized  and  teflon  coated. 

3.1.6  Valves 

Cartridge  type  brass  check  valves  in  stainless  steel  housing. 

3.1.7  Intercoolers 

Stainless  steel  tubing  intercoolers. 

3.2  Design  Philosophy 

The  purpose  of  building  the  prototype  was  to  obtain  in  as  short  a 
time  as  possible,  preliminary  results  on  the  performance  of  all  the 
key  components  under  normal  operating  circumstances.  Off-the-shelf 
components  were  used  where  ever  possible. 

The  design  of  the  piston  seals  called  for  three  sealing  rings  per 
piston,  however,  only  single  seals  were  installed  to  allow 
measurements  of  seal  wear  at  an  accelerated  rate. 

Engineering  drawings  have  been  completed  of  the  cylinders,  piston 
and  valve  assemblies.  These  have  been  fabricated  and  assembled  and 
put  through  preliminary  tests. 

3.3  Performance 

The  prototype  compressor  met  the  design  goal  delivering  2.0  1/min 
(STP)  at  2000  psi  outlet  pressure.  At  zero  flow  the  outlet  pressure 
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Compressor  Schematic 


exceeded  3200  psi ,  limited  by  the  hydraulic  pressure  relief  valve 
setting. 


Approximately  100  hours  of  testing  of  the  compressor  was  completed 
on  dry  air  with  a  moisture  content  below  lppm.  The  compressor  was 
then  operated  in  a  closed  cycle  mode  with  a  MMR  microminiature 
refrigerator  and  gas  drier  in  the  circuit.  The  test  was  interrupted 
periodically  to  inspect  the  check  valves,  sealing  surfaces  and 
hydraulic  surfaces.  A  total  of  approximately  370  hours  of  operation 
has  been  logged. 

3.4  Results  and  Discussion 

The  only  major  problem  encountered  was  leakage  of  hydraulic  fluid 
past  the  rod  seals  on  the  hydraulic  cylinder.  The  cylinder  used, 
was  a  low  cost,  low  pressure  cylinder  which  was  available 
off-the-shelf.  It  was  operated  at  800  psi  which  is  close  to  its 
maximum  rating.  Discussions  with  the  manufacturer  indicates  that 
use  of  a  higher  rated  cylinder  and  one  with  double  rod  seals  will 
eliminate  this  problem. 

One  failure  did  occur  with  a  high  pressure  seal.  However, 
examination  of  the  seal  indicates  that  it  was  probably  nicked  during 
assembly  and  failed  to  seat  upon  being  pressurized.  The  seal  had 
clearly  not  worn  significantly  but  had  folded  over  and  broken. 


The  prototype  has  proved  the  general  validity  of  the  approach.  The 
neccesary  changes  are  being  made  in  the  choice  of  hydraulic  cylinder 
and  some  redesign  is  being  made  of  the  end  plates  to  strengthen  the 
support  of  the  valve  assemblies.  A  second  prototype  will  then  be 
assembled  and  operated  over  an  extended  period  to  check  on  the  wear 
of  tne  sealing  rings  and  check  valves  before  proceeding  with  the 
addition  of  the  second  and  third  gas  circuits. 

4.0  Objectives  of  the  Next  Reporting  Period 

During  the  next  reporting  period  it  is  expected  that  the  following 
tasks  will  be  completed: 

Experimental 

1.  The  first  five  layers  of  the  refrigerators  will  be  fabricated. 
Figure  4.0-1  '•hows  an  example  of  the  "two  by  two"  step  and 
repeat  patt-  s  that  are  required  for  the  photo-resist  etching 
process.  The  manufacture  of  the  glass  components  for  the  first 
five  layers  will  begin  with  this  step. 


Each  stage  (argon  and  hydrogen)  will  be  assembled  independently 
and  flow  tested. 

The  argon  stage  will  be  tested  for  full  cryogenic  performance 
with  high  pressure  gas. 
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4.  Once  the  argon  stage  is  found  to  operate  properly,  the  hydrogen 
stage  will  be  attached  and  the  assembly  will  be  operated  with 
high  pressure  gas. 

5.  A  second  single  gas  prototype  compressor  will  be  completed  and 
test  on  it  will  be  initiated.  This  will  then  be  used  for  tests 
of  each  of  three  gases  prior  to  the  fabrication  of  the  three  gas 
unit. 

Theoretical 

1.  A  full  layout  of  the  helium  stage  will  be  completed. 

2.  Photo-masks  will  be  made  to  prepare  for  fabrication  of  this 
stage. 


5.0  Glossary 

The  variables  used  in  the  calculations  are  defined  as  follows 


q  = 
q  x 

m  = 
m  = 
H  = 
K  = 
t  = 
r  = 
R  = 
L  = 
T  = 
C  = 
AH  = 
^T(t)  = 

aTc  * 
h  = 

l  = 
w(0  = 

V  = 


u  = 

V  = 


p  = 


C  (T) 
Ptc 


d  = 


o  = 


e 


thermal  mass,  joules 
power  watts 
mass  grams 

mass  flowrate,  gm/sec 
enthalpy,  joules/gr 
thermal  conductivity,  watts/cm  K 
thickness,  cr. 
radius,  cm 

electrical  resistance,  ohms 

-8  2 

The  Lorentz  No.  =  2.45  x  10  (volts/K) 

temperature,  K 

refrigeration  capacity,  vatts 
enthalpy  difference  joules/err 
time  dependent  temperature  difference  (K) 
initial  temperature  difference  (K) 
thermal  diffusivity  (cm  /sec) 
characteristic  length  (cm) 

width  of  refrigerator  disk  as  function  of  length 
volumetric  flowrate  at  STP,  liter/min 
Jcule-Thomson  Coefficient  (h/atm) 
volume,  cm 

3 

density,  gm/cm 

temperature  dependent  specific  heat,  (joules/gm  K) 
time,  min. 
diameter,  cm 

-12  2 

Stefan  Boltzmann  Constant,  5.67  x  10  watts/cm  K 
emissi vity 

-1 

electrical  conductivity,  (ohr  cm) 
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R 

o 

Mk  = 

P  * 
o 

a  = 
h  = 


gas  constant  8. 2  x  10  ,  (liter  atm/g  mole  K) 
molecular  weight,  (gm/g  mole) 
exhaust  pressure,  (atm) 

refrigerator  flow  coefficinet,  (liter/atm  min) 
heat  transfer  coefficient,  watts/cm  K 
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